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The digitized readout (DN = data number) of the low-voltage signal (MV) which 
is proportional to the potential (V) on the high-voltage modulator grids, is a part of the 
routine calibration sequence on the Voyager Plasma Science (PLS) experiment. 
Theoretically, in producing the modulator voltage, the voltage is stepped through one 
decade from 60 to 600 volts and then through the same steps with a multiplication fac- 
tor of 10 to sweep from 600 to 6000 volts. To this voltage is added an offset of -50 
volts. 

The voltage monitors measure a voltage between 0.1 and 10 volts (a constant 
fraction of the voltage swept) which is then input into the fast A/D converter. This 
voltage is converted to a binary number from which the highest order bit is discarded. 
Thus, for MV sweeping from 0.1 to 1.0, a number between 0 and 255 is returned; for 
MV sweeping from 1.0 to 10 volts, again, a number between 0 and 255 is returned. 
To know which decade is being read out, one must either rely on the “proper 
sequence” or an “educated guess.” 

For more information on MVM Interpretation, see Voyager Memorandum #161 
by R. L. McNutt, Jr. (attached). 


2. An Outline of the PLS Calibration Data Analysis 

The PLS Modulator Calibration (MVM) Data Analysis was undertaken in order to 
check the correctness of the fast A/D converter formulas that connect low-voltage 
monitor signals (MV) with digital outputs (DN), to determine the proportionality con- 
stants between the actual modulator grid potential (V) and the monitor voltage (MV), 
and to establish an algorithm to link the digitized readouts (DN) with the actual grid 
potential (V). The data used for the DN-MV analysis were from the calibration tests 
for the PLS-Prototype run at MIT CSR (included as Appendix A of Voyager Memo 
#161). The MV-V analysis used the results of the power supply tests, made by Matrix 
Research and Development Corp. for both PLS instruments to be flown on the Voy- 
ager 1 and Voyager 2 spacecraft (in the further text referred to as “test results”). The 
Matrix data was obtained between 11/11/76 and 11/15/76 for “Flite 1,” and between 
12/21/76- and 12/27/76 for “Flite 2.” Drawings containing this data are labeled 
“MJ.S. Power Supply” (there are no drawing numbers). These are stored in the filing 
cabinet next to the outer wall, second drawer from the top in N5 2-367. Note that there 
is some potential confusion about which unit is which: if Flite 1 is SN002, then it is 
the instrument actually flown on Voyager 1. SN002 was put on Voyager 1 and SN003 
which was supposed to have gone on Voyager 2 was not flown. When these data were 
taken, Flite 2 probably referred to SN003 which has remained at M.I.T. as the flight 
spare. SN001 was actually flown on Voyager 2. 
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All calibration tests and power supply tests were entered into a sequence of files. 
The method that was used in both analyses, i.e., DN-MV and MV-V, is an improved 
version of a least squares fit algorithm, implemented in several FORTRAN 77 pro- 
grams (see, for example section 14.2 of Numerical Recipes: The Art of Scientific Com- 
puting, Press, W. H., B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, Cambridge 
University Press, New York, 1986). 

3. Description of the Files 

All of the files related to the investigation can presently be found on the directory 
/usr4/ojn. The file suffixes have the following meanings: 

.dat data files 

.mon mongo source code files 

.f FORTRAN 77 source code files 

The DN-MV test data and mongo files have filenames of the form: 
p[ mode ][ mvr # ].[ suffix ] 

The MV-V test data and mongo files have filenames of the form: 
f[ flite # ]-[ temperature ]-[ mode ].[ suffix ] 


where: 


mvr # = 1 if 0.1 £ MV £ 1.0 
mvr # = 2 if 1.0 £ MV <, 10 

flite # = 1 for Voyager 1 PLS Instrument 
flite # = 2 for Voyager 2 PLS Instrument 

temperature = 5 for Temp = -5°C 
temperature = 10 for Temp = 10°C 
temperature = r for Temp = room-temp 
temperature = 30 for Temp = 30°C 
temperature = 40 for Temp = 40°C 

(The temperature is that to which the instrument was exposed during the test.) 

mode = 1 for the L-mode 
mode = m for the M-mode 
mode = el for the El-mode 
mode = e2 for the E2-mode 
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suffix = dat if the file is a data file 

suffix = mon if the file is a mongo source code 


Note that there are no files of the form f2-30-[ mode ].[ suffix ] due to lack of 
this information. 

The files of the form: 
tl[ mvr # ].f and tl[ mvr # ].dat 

are FORTRAN 77 source code and output data files produced by their corresponding 
source code programs which make/contain the data from the theoretical fast A/D for- 
mulas for MV-DN conversion. 

The FORTRAN 77 source code files that actually perform the analysis are: 

dns-code.f DN-MV analysis 

mvs-code.f MV-V analysis 

mvs-tc-code.f MV-V analysis combining the data files 

that differ only by temperature 


These source code files use the following driver files: 

dns-driver.dat 

mvs-driver.dat 

mvs-tc-driver.dat 


The analysis generates the following files: 


dns-output.dat 

vs-output.dat 

mvs-tc-out.dat 


DN-MV analysis results 

MV-V analysis results 

MV-V temperature-combined analysis results 
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4. Data File Organization 

All data files are organized as matrices of numbers and/or strings. Their row 
structure is as follows: 

p[ mode ][ mvr # ].dat and tl[ mvr # ].dat: 
dn mv 

where: dn are the data numbers 

mv are the corresponding monitor voltages 

f[ flight # ]-[ temperature ]-[ mode ].dat: 

channel# hvdcout mvu mvl vl vu 
where : channel# is the channel number 

hvdcout is the dc high voltage output 
mvl/mvu are the monitor low voltages at the 
lower/upper channel edge 
vl/vu are the actual grid potentials corresponding 
to the mvl/mvu 

dns-driver.dat: 

filename mvr[ mvr # ] lines# 
where : filename is the name of a file. 

mvr[ mvr # ] is a character string that can be either 

mvrl or mvr2 depending on what range the value 
of MV belongs to. 

lines# is the number of lines in the file. 
mvs-driver.dat: 

filename startline- 1 endline- 1 startline- 2 endline-2 

where : startline- [ chr # ]/endline-[ chr # ] are the 
starting/ending line numbers for parts of the 
f[ ]-[ ]-[ ].dat data files that correspond to different 
channel ranges as follows: 

chr # = 1 if 0 £ channel# £ 7 and mode = L or E2 

0 £ channel# £ 63 and mode = M 

all channels and mode = El 

chr # = 2 if 8 < channel# £ 15 and mode = L or E2 
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64 £ channel# £ 127 and mode = M 

mvs-tc-driver.dat: (Note: this file is a single-column file.) 

mfilename#l 
filename# 1.1 

filename# 1.2 
* 

mfilename#2 

filename#2.1 

where : mfilename#[ integerl ] is a character string of the form: 
f[ flite # ]-[ mode ].dat 

which is found in the mvs-tc-outpuLdat file and 
labels the analysis results obtained by combining 
groups of f[ ]-[ ]-[ ].dat files that differ only by 
temperature. 

filename#[ integerl ].[ integer2 ] are the filenames of the 
f[ ]-[ ]-[ ].dat data files that were combined. 

dns-output.dat: 

filename mvr[ mvr # ] cl c2 sqd pairs# 

where : filename is the name of a file that labels the information 

obtained by its analysis. 

mvr[ mvr # ] is as explained under dns-driver.dat. 

Ci, c 2 are constants in the presumed function 
that fits the data. 

sqd is the square deviation of the function, 
pairs# is the number of analyzed data pairs. 

mvs-output.dat: 

filename chrf chr # ] mvx-vx a b sqd siga sigb cov pairs# 

where : chr[ chr # ] is a character string that can be either 

chrl or chr2. chr # has a meaning as explained 
under mvs-driver.dat 

mvx-vx is a character string that can be either 

mvl-vl or mvu-vu. It specifies what monitor 
voltages were analyzed - those at the upper or 



i 


- 6 - 


at the lower channel edge. 
a, b are constants in the presumed function that 
fits the data. 

sqd, siga (<jJ, sigh (a i) )are, 

respectively, the square deviation 
of the function, and the deviations of its 
coefficients a and b. 
cov (Gab) is the covariance of a and b. 

mvs-tc-output.dat 

mfilename chr[ chr # ] mvx-vx a b sqd siga sigb cov pairs# 

where : mfilename is the character string that labels the analysis 

data, and contains information what files were 
combined for the analysis (see mvs-tc-driver). 


5. Mongo Files 

All files that have a suffix .mon are mongo source code files. There are two types 
of such files: p[ ][ ].mon and ft ]-[ ]-[ ].mon files. 

Every mongo file is used to create a hardcopy plot of the data file associated with 
the mongo file. A plot is obtained by typing: 

mongo [ mongo-filename ] 

p[ mode ][ mvr # ].mon plots include both the p[ mode ][ mvr # ].dat data (square 
dots) and the tl[ mvr # ].dat data (broken line). f[ flite # ]-[ temperature ]-[ mode 
].mon files produce four plots out of the associated f[ ]-[ ]-[ ].dat data for different 
combinations of the channel range and choice of mvl-vl or mvu-vu pairs of data 
(square dots). 


6. Fortran Files 

Files that end with .f are FORTRAN 77 source code files. There are five such 
files: tll.f, tl2.f, dns-code.f, mvs-code.f, and mvs-tc-code.f. The sequence of UNIX 
commands: 

f77 -o execfile filename.f 
execfile 



k 


- 7 - 

i 

I 

performs a compilation of the source code and an execution of execfile, the gen- 
erated object file. 

tll.f & tl2.f: 


These programs generate the theoretical DN-MV pairs assuming the following 
exponential functions: 


\fy — e -(255-DN)tft 


for tll.f i.e. mvr# = 1 

MV = 10 e - < 255-£W ^ / ' c 


for tl2.f, i.e., mvr# = 2 


( 1 ) 

( 2 ) 


where t = 4.34 (is, x = 482.55 ps and 0 £ DN <, 255 (From A. Mavretic’s lab note- 
book labeled book 2, #254, July 1974). Note that = 10 1 ' 256 0167 for these values. 

dns-code.f : 

This program does a least squares fit by analyzing the p[ ][ ].dat data files, and 
presuming that MV and V are related by the following exponential relation: 

MV = ci e ~ { ' 25S ~ DN) Cl (3) 

where c x and c 2 are constants that are determined by the fitting procedure. 

The program produces the file dns-output.dat. The algorithm used is a 
modification of the general least squares algorithm that makes the round-off errors as 
small as possible. 

mvs-code.f : 

This program uses the same algorithm as the previous one. It presumes that MV 
and V are connected with the following relation: 

V = a MV + b (4) 

where a and b are constants that are determined by the fit. 

The program analyses the f[ ]-[ ]-[ ].dat data files, dividing each file into four 
parts using the driver file. Every part is characterized by a unique channel range and 
mv range (mvr). The analysis, i.e., the least squares fit is carried out independently for 
each part. The file mvs-output.dat contains the results from the analysis. 

The formulae used to calculate the deviations of the obtained coefficients are: 


i 
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1 1/2 



(5) 

a b = [v a 2 / a] 

(6) 

dab = “X MV i a 2 / A 
i 

(7) 


where 


and 


A*N£MV?- 


£ MV, 

i 




° 2s 1^2 Sqd (8) 

with N the number of MV-V pairs, 
mvs-tc-code.f : 

This program employs the same least squares fit algorithm as the previous two. It 
presumes the same relation between MV and V, and performs the same analysis as 
mvs-code.f does, except that it internally combines the f[ ]-[ ]-[ ].dat data files that 
differ only by temperature. The results of the analysis are contained in the file mvs-tc- 
out.dat. 


7. DN-MV Analysis 

The primary goal of this analysis was to find out how much the theoretical DN- 
MV formulas are off from the true values. The only information available is from the 
calibration tests performed on the PLS-Prototype. This information was severely lim- 
ited - only one test, at unknown temperature, per channel in each mode. Also, there is 
no information on the measurement uncertainties. Hence, the analysis results should 
be used with some caution, consideration. 

The mongo plots of DN-MV data files (DN vs. log MV) show that DN and MV 
are related through an exponential relation. The DN-log(MV) pairs, as can be seen 
from the plots, lie on a straight line which, however, has a slightly different slope and 
y-axis intercept than the theoretical line. This suggested that there might be some non- 
linearity in the A/D conversion, if one assumes that DN’s and MV’s were measured at 
the same time. In addition, the slope and intercept shifts are bigger for mvrl (0.1 
< MV <. 1.0) than the corresponding shifts for mvr2. This implies that the multiplica- 
tion factor in the DN-MV formula for this range is not exactly 10. 

The investigation of the plots sugested that the presumed function for the least 
square fit be: 
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MV = ci e H25S DN) Ci (9) 

The program actually uses a linear least squares fit algorithm. If we take In of the 
both sides of the equation, we obtain: 

In MV=(DN-255) c 2 + In (10) 

which can be considered as a linear function y = ax + b setting: x = DN- 255; 
y = In MV ; a = c 2 ; b = In 

Before the data files were used for the analysis, all DN-MV pairs which were 
significandy off the imagined line that connects the other points on the mongo plots 
were deleted. They obviously represent a subjective reading, or even a writing, error. 

The analysis results show that the presumed function yields a good approxima- 
tion, since the square deviation is less than 10 -6 in most cases. The multiplication fac- 
tor cj ranges from: 

1.02831 -► 1.02864 instead of 1 for mvrl (2.831% -> 2.864%) 

9.91158 — > 9.92186 instead of 10 for mvr2 (-0.8842% -0.7814%) 

The range deviations for c t in percentages are: 

mvrl: ± 0.016% 
mvr2: ± 0.052% 

is surprising that the c 1 (mvr2)/c 1 (mvr2) ratio is less than 10 ; it ranges from 9.63880 
-» 9.64560. This implies that PLS circuitry constants responsible for the multiplication 
factor might be off by -3.5% from their nominal value(s). 

The constant in the exponent, c 2 = t / x, ranges from : 

0.888580 x 1(T 2 -» 0.889623 x 10 ~ 2 instead of 0.899388 x 10“ 2 (1.201% 
1.085%). 

The range deviation for c 2 is ± 0.058%. 

If we consider t (the period of the fast A/D converter) as constant, then the time 
constant x might be off (bigger) 1.085% — » 1.201% from the nominal value of 482.55 
p.s, or, on an average, by - 1.15%. 
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8. MV-V Analysis 

The goal of this analysis was to determine the proportionality constants between 
the monitor low voltage and high voltage (i.e., the actual grid potential). The low vol- 
tage is (nominally) proportional to the high voltage applied to the modulator grids (at 
least in the linear regime of the amplifiers). The theoretical value of this factor is not 
known, although it can, in principle, be derived from the circuit diagrams. 

The analysis uses the results of the power supply test, conducted by Matrix 
Research and Development for two of the PLS instruments (refer to discussion in sec- 
tion 2). The tests were performed for several temperatures. Although all modes were 
incorporated in the measurements, the files were highly incomplete and difficult to read 
(handwritten). As is the case with the DN-MV analysis, there is no information on the 
measurement uncertainties. 

All of the plots of MV-V data files clearly show that MV and V are linearly 
dependent. Since no information was readily available for the theoretical factor of pro- 
portionality, the plots do not include any theoretical graphs. A simple linear function 
was presumed for the least squares fit: 

V = a MV + b (11) 


Similar to the DN-MV analysis, all MV-V points with unusually large deviations 
from this function were deleted from the files. 

Throughout the analysis, every file is divided into four parts (with an exception of 
mode El files which are divided into two parts with respect to mvr). Each part has a 
unique channel range (chrl or chr2), and monitor monitor voltage range (mvrl or 
mvr2). These parts are processed independently. 

The analysis results (file mv-outputdat) are somewhat surprising. First, the square 
deviation ranges from roughly 10 -5 to over 350 (considering least squares fits of more 
than 2 data pairs). This might be evidence of harsh measurement errors during the test 
or strong non-linear influences in the PLS circuit The random distribution of the 
square deviation values shows that the errors or influences are not dependent on factors 
such as: temperature, mode, voltage or channel range. The deviations obtained for a 
and b, as well as the covariance of a and b, show a much more stable picture; their 
median values are, roughly, as follows: 

c a -0.1 
c b - 0.1 
Gab ~ “IQ -2 
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It is interesting that the values obtained for a and b show an obvious dependence on 
various factors, although, ideally, they should not. If one looks at the values for the 
slope a, it can be seen that in 90% of the cases: 

a(mvl-vl) > a(mvu-vu) 

There are only eight cases where this relation does not hold true; since these cases 
appear to be totally random with respect to any factors that might have influenced 
them (temperature, specific mode, etc.), the observed property may serve as an evi- 
dence that some circuit constants change when switching from the lower to the upper 
channel edge. 

The ratio a(chr2)/a(chrl) randomly oscillates around the value of 10, which is the 
nominal value for the multiplication voltage factor when the instrument switches from 
the first channel range to the second channel range. The stability of this multiplication 
factor shows that the circuit constants for the appropriate part of the circuit are close to 
their nominal values. 

This situation, however, changes in considering the dependence of the slope of a 
on the mode. The values of a for various modes vary roughly as follows: 

1-mode a = 60.5; 605.1 

m-mode a = 60.3; 603.5 

el-mode .... a = 67.2 
e2-mode .... a = 68.1; 672.1 

The slight difference between a(l-mode) and a(m-mode) can be noticed in many cases, 
especially considering values of a which correspond to a smaller square deviation and 
larger number of analyzed MV-V pairs. On die average, a(l-mode) is 0.33% larger 
than a(m-mode). Following the same criterion, one notices that a(el-mode) and a(e2- 
mode) differ more drastically from a(l-mode) and a(m-mode). The percentage increase 
from a(m-mode) turns out to be, roughly: 

a(el-mode) 11.45% 

a(e2-mode) 12.93% 

This suggests that some circuit constants in both E modes differ from their nominal 
values. 

The last consideration for the slope values is whether or not they are temperature 
dependent It is important to note, once again, that the analysis was limited to five 
different temperatures. The values of a do not show any dependence on temperature; 
they tend to oscillate randomly as the temperature increases. This observation is con- 
sistent with the theoretical expectations. 
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The values of b, which represent the y-axis intercept or the voltage offset are 
much more randomly distributed than are the slope values. In fact, b values vary so 
much that no dependence property can be uniquely established. 

The values of b vary from -40 to +40; however, most of them are between -5 to 
0. In roughly 90% of cases the values are negative. Theoretically, the values of b 
should be close to 0. Typical values for b are (roughly): 

-0.4 ... for chrl 

-1.5 ... for chr2 


9. MV-V Temperature-combined Analysis 

Since the slope values did not show any dependence on temperature, all families 
of files that differed only by temperature were collapsed, and the data from the 
corresponding parts of these files combined for the least squares fit. The result of this 
analysis is the file mvs-tc-outdat. First, it can be noticed that the lower edge of the 
square deviation range has moved toward larger values (from -10 -5 to ~10 -3 ). The 
square deviations for a and b, and the covariances are, for the most part, in the same 
range as in the previous analysis. 

The relation: 

a(mvl-vl) > a (mvu-vu) 

still holds true for 80% of cases (there are three cases when it does not). 

The ratio a(chr2)/a(chrl) oscillates around 10 (the theoretical value) as in the 
MV-V analysis. 

This analysis also confirms that there is a definite dependence of a on the mode. 
Some spot checks suggest the following percentage differences (with respect to a(m- 
mode)): 

a(l-mode) 0.13% 

a(e 1-mode) .... 11.75% 
a(e2-mode) .... 12.12% 

(The median value for a(m-mode) is 60.13.) 

The values of b are, again, very randomly distributed, so that no firm property or 
dependence can be established. For example, b me = -4.631 volts. However, if the first 
four large values of b are not considered (since they range from -20 to -30), the new, 
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restricted average value of b becomes: b me ra = -1.153 volts. 

This example clearly shows that it is very difficult to derive any major conclusion 
about b from the analysis results obtained. 


10. Conclusions 

The analysis results are surprising in that the derived conversion constants deviate 
by fairly significant amounts from their nominal values. However, it must be kept in 
mind that the test results which were used for analysis may be very imprecise. Even if 
we assume that the test result errors are very large, they do not appear to be capable to 
account for all discrepancies between the theoretical expectations and the results of the 
analysis. Measurements with the flight spare instrument appear to be the only means 
of investigating these effects further. 

It is very clear that is impossible to create one simple algorithm that for given 
DN will return V - the actual grid potential. The MV-V slope dependence on the chan- 
nel range is what was expected. However, the MV-V slope dependence on whether 
we have the top or the bottom of the channel, and a dependence of the MV-V slope on 
the mode, precludes a unique algorithm for all cases. 

In order to convert a given DN into the V, we need the following auxiliary infor- 
mation: 

1 - the mode 

2 - the decade the modulator is sweeping through (mvr#) 

3 - the channel range (chr#) 

4 - the channel edge the DN corresponds to (mvl-vl or mvu-vu) 

Given the stated information (for 2 and 3 we must rely on the “proper sequence” or 
an “educated guess”), one can pick up the corresponding constants from the files 
dns-output.dat and mvs-te-out.dat, and substitute them in the general formula: 

V = a c x (255 " DW) + b (12) 

to obtain the actual grid potential V. The error in the value obtained will be as large as 
the errors in the test results the DN-MV-V analysis is based on. 

For example: Flite# = 2 ; Mode = M ; mvr# = 2 ; chr# = 1 ; channel edge = 
lower (mvl-vl). The values of the constants in the general DN-V formula, are: 

c x = 9.91168 c 2 = 0.888609 x 10“ 2 
a = 60. 13527 b = 0.71937 x 10" 1 

If the value returned from the PLS instrument is DN = 73, then the grid potential is 
118.3 volts. 
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It is important to note that the general DN-V formula only yields magnitudes; the 
actual potentials in the El and E2 modes will be negative. 


11. Description of Tables 

Tables 1,2, and 3 are, respectively, hardcopies of the files dns-output.dat, mvs- 
output.dat, and mvs-tc-out.dat. 


12. Description of Graphs 

The first set of graphs contains hardcopy plots of p[ ][ ].mon mongo files. The 
second set of graphs contains hardcopy plots of f[ ]-[ ]-[ ].mon mongo files (see sec- 
tion 5). 


13. Description of Source Code used in the Analysis 

Hardcopy listings (with comments) of the files dns-code.f, mvs-code.f, and mvs- 
tc-code.f are included. 

14. Voyager Memorandum #161 
Voyager Memorandum #161 is appended. 
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DNS-CODE.F analyses the PLS Prototype data. 
Author Ognen Nastov , MIT undergraduate, Nov 1987 


leastsq does a least square approximation 

and minimizes the round-of errors, 
n » # of x-y pairs 
A - (1 / sum(t”2)) * sum (t’y) 

B » (sum(y) - A’sum(x)) / n 
where t(i) - x(i) - sum(x) / n 
subroutine leastsq (filename, depth, 
a, b, sqd) 
integer depth, dn 
real a, b, mv, x, y, sumx, sumy, 

t(120), u(120), v(120), sumtsq, sumty, sqd 
character*8 filename 
sumx « 0 
sumy - 0 
sumtsq - 0 
sumty - 0 

open (1, file-filename, form- ’formatted’) 
rewind (1) 

do 1010 counter - 1, depth 
read(l, *) dn, mv 
x - 255 - dn 
y - log (mv) 
sumx - sumx + x 
sumy - sumy + y 
u(int(counter)) - x 
v(int(counter)) - y 
continue 

do 1020 i - 1, depth 

t(i) - u(i) - sumx / depth 
sumtsq - sumtsq + t(i)*t(i) 
sumty - sumty + t(i)*v(i) 
continue 
a - sumty / sumtsq 
b - (sumy - a’sumx) / depth 
call sqdev (a, b, depth, u, v, sqd) 
close (1) 
return 
end 


c Sqdev calculates the square deviation 

c by definition. 

subroutine sqdev (a, b, n, u, v, sqd) 
integer n 

real a, b, u(120), v(120), sqd 
sqd-0 

do 2000 i - l,n 

sqd - sqd + (a*u(i) + b - v(i))”2 
2000 continue 

sqd - sqd/ n 
return 
end 


c The MAIN program body, 

c It produces one line of output for each input file, 

c The program uses a special driver file 

c dns-driver.dat. 

integer pairs 
real a, b, cl, c2 
character’s filename 
character’4 mvr 

open (2, file- ’dns-driver.dat’ , form- ’formatted’) 
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open (3, file- ’dns-outputdat’ , form- ’formatted’) 
rewind (2) 

print *, 'Creating file: ’, ’dns-output.dat’, 
c 

50 read (2, *) filename, mvr, pairs 

print *, ’Processing file: ’, filename, 
call leastsq (filename, pairs, a, b, sqd) 
cl = exp (b) 
c2 = a 

write (3, fmt-200) filename, mvr, cl, c2, sqd, pairs 
100 if (filename .ne. ’pe22.dat’) go to 50 
close (2) 
close (3) 

print *, ’File: dns-outputdat created.’ 

200 format (a8, tr2, a4, tr2, gl2.6, tr2, gl2.6, tr2, 

+ gll.6, tr2, i3) 


end 
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c MVS-CODE.F analyses the PLS calibration data for both flites. 

c 

c Author. Ognen Nastov, MIT undergraduate, Nov 1987 

c leastsq does a least square approximation 

c | sum(x**2)A + sum(x)B - sum(xy) 

c j sum(x)A + nB - sum(y) 

c where n - # of x-y pain 

subroutine leastsq (filename, voltrange, startline, endline, 

+ a, b, sqd, depth) 

integer channel, startline, endline, skip, depth 
real a, b, mvu, mvl, vl, vu, hvdcout, x, y, sumxsq, sumx, sumxy, 
+ sumy, sumysq, sqd 

character"^ voltrange 
character* 12 filename 
sumxsq - 0 
sumx - 0 
sumxy - 0 
sumy - 0 
sumysq - 0 

open (1, file-filename, form- ’formatted’) 

rewind (1) 

skip - startline • 1 

do 500 counter - 1, skip 

read (1,*) channel, hvdcout, mvu, mvl, vl, vu 
500 continue 

depth - endline • startline + 1 
do 600 counter - 1, depth 

read (1, *) channel, hvdcout, mvu, mvl, vl, vu 
if (voltrange .eq. ’lower’) then 
x » mvl 
y - vl 
else 

x - mvu 
y - vu 

endif 

sumxsq - sumxsq + x*x 
sumx - sumx + x 
sumxy - sumxy + x*y 
sumy - sumy + y 
sumysq - sumysq + y*y 
600 continue 

call solvesys (sumxsq, sumx, sumxy, sumx, depth, sumy, 

+ a,b) 

call sqdev (a, b, depth, sumx, sumy, sumxy, 

+ sumxsq, sumysq, sqd) 

close (1) 
return 
end 


c 

c 

c 

c 

c 

c 


+ 

+ 

+ 


lsq2 does a least square approximation 

and also minimizes round-of errors 
n - # of x-y pairs 
A - (1 / sum(t**2)) * sum (t*y) 

B - (sum(y) - A*sum(x)) / n 
where t(i) = x(i) - sum(x) / n 
subroutine lsq2(filename, voltrange, startline, endline, 
a, b, sqd, sqda, sqdb, cov, depth) 
integer channel, startline, endline, skip, depth 
real a, b, mvu, mvl, vl, vu, hvdcout, x, y, sumx, sumy, 
t(20), u(20), v(20), sumtsq, sumty, sqd, 
sqda, sqdb, cov, del, sumxsq 
character*6 voltrange 
character* 12 filename 
sumx - 0 
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sumy - 0 
sumtsq - 0 
sumty - 0 
sumxsq « 0 

open (1, file-filename, form- ’formatted’) 

rewind (1) 

skip - startline • 1 

do 1000 counter - 1, skip 

read (1, *) channel, hvdcout, mvu, mvl, vl, vu 
1000 continue 

depth - endline - startline + 1 
do 1010 counter - 1, depth 

read (1, *) channel, hvdcout, mvu, mvl, vl, vu 
if (voltrange .eq. ’lower’) then 
x - mvl 
y = vl 
else 

x- mvu 
y- vu 

endif 

sumx - sumx + x 
sumy - sumy + y 
sumxsq - sumxsq + x*x 
u(int(counter)) - x 
v(int(counter)) - y 
1010 continue 

do 1020 i - 1, depth 

t(i) - u(i) - sumx / depth 
sumtsq - sumtsq + t(i)*t(i) 
sumty - sumty + t(i)*v(i) 

1020 continue 

a - sumty / sumtsq 
b - (sumy - a*sumx) / depth 
call sqd2 (a, b, depth, u, v, sqd) 
close (1) 

del - depth * sumxsq - (sumx)**2 

sqdb - sqrt (sumxsq / del * sqd * depth / (depth - 2)) 

sqda » sqrt (depth / del * sqd * depth / (depth - 2)) 

cov - - sumx / del * sqd *depth / (depth - 2) 

return 

end 

c sqd2 calculates die square deviation by definition 

c sqd - (1 / n) * sum(Ax + B - y)**2 

subroutine sqd2 (a, b, n, u, v, sqd) 
integer n 

real a, b, u(20), v(20), sqd 
sqd»0 

do 2000 i- 1, n 

sqd - sqd + (a*u(i) + b - v(i))**2 
2000 continue 

sqd - sqd / n 
return 
end 

c solvesys solves 2x2 linear systems using determinants 

c al*x+bl*y=cl 

c a2*x+b2*y-c2 

subroutine solvesys (al, bl, cl, a2, b2, c2, x, y) 
integer b2 

real al, bl, cl, a2, c2, x, y, det, deu, dety 
det - al*b2 -a2*bl 
detx - cl*b2 - c2*bl 
dety » al*c2 - a2*cl 
x - detx/det 
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y = dety/det 

return 

end 

c sqdev calculates the square deviation expanding 

c the sum: 

c sqd- (1 /n) * sum(Ax + B - y)**2 

subroutine sqdev (a, b, n, sumx, sumy, sumxy, 

+ sumxsq, sumysq, sqd) 

integer n 

real a, b, sumx, sumy, sumxy, sumxsq, sumysq, sqd 
sqd = a*a*sumxsq + n*b*b + sumysq + 2*a*b*sumx - 
+ 2*a*sumxy - 2*b*sumy 

sqd = abs (sqd / n) 
return 
end 

c The MAIN program body, 

c It produces four lines of output for each input file, 

c The program makes use of the lsq2 and sqd2 subroutines, 

c It uses the file mvs-driver.dat to get the info 

c essential to analyse only the appropriate chunks 

c of data. 

integer stll, edll, stl2, edl2, pairsl, pairsu 
real al, bl, sqdl, sqdal, sqdbl, covl, 

+ au, bu, sqdu, sqdau, sqdbu, covu 

character*6 voltrange 
character* 12 filename 

open (2, file- ’mvs-driver.dat’, form- ' formatted’ ) 
open (3, file-’mvs-outputdat’, form- ’ formatted’ ) 
rewind (2) 

print *, ’Creating file: ’, ’mvs-outputdat’, ’...’ 
c analyze the first channel range... 

50 read (2, *) filename, stll, edll, stl2, edl2 
print *, ’Processing file: ’, filename, ’...’ 
voltrange - ’lower’ 

call lsq2 (filename, voltrange, stll, edll, 

+ al, bl, sqdl, sqdal, sqdbl, covl, pairsl) 

voltrange - ’upper’ 

call lsq2 (filename, voltrange, stll, edll, 

+ au, bu, sqdu, sqdau, sqdbu, covu, pairsu) 

write (3, fmt-200) filename, ’chrl’, ’mvl-vl’, al, bl, 

+ sqdl, sqdal, sqdbl, covl, pairsl 

write (3, fmt-200) filename, ’chrl’, ’mvu-vu’, au, bu, 

+ sqdu, sqdau, sqdbu, covu, pairsu 

c analyze the second channel range... 
if (stl2 .eq. 0) go to 100 
voltrange - ’lower’ 

call lsq2 (filename, voltrange, st!2, ed!2, 

+ al, bl, sqdl, sqdal, sqdbl, covl, pairsl) 

voltrange - ’upper’ 

call lsq2 (filename, voltrange, stl2, edl2, 

+ au, bu, sqdu, sqdau, sqdbu, covu, pairsu) 

write (3, fmt-200) filename, ’chr2\ ’mvl-vl’, al, bl, 

+ sqdl, sqdal, sqdbl, covl, pairsl 

write (3, fmt-200) filename, ’chr2’, ’mvu-vu’, au, bu, 

-i- sqdu, sqdau, sqdbu, covu, pairsu 

100 if (filename .ne. ’f2-40-e2.dat’) go to 50 
close (2) 
close (3) 

print *, ’File: mvs-output.dat created.’ 

200 format (al2, 02, a4, tr2, a6, tr2, f9.5, tr2, gl2.6, tr2, 

+ gl 1.6, tr2, gl2.6, tr2, gl2.6, tr2, gl2.6, tr2, i2) 


end 
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c MVS-TC-CODEJ 

c Author Ognen Nastov, MIT undergraduate, Nov 1987 

c lsq2 does a least square approximation 

c and also minimizes round-of errors 

c n « # of x-y pairs 

c A - ( 1 / sum(t**2)) * sum (t*y) 

c B • (sum(y) - A*sum(x)) / n 

c where t(i) - x(i) - sum(x) / n 

subroutine lsq2(filenames, filenum, voltrange, startlines, 
+ endlines, a, b, sqd, sqda, sqdb, cov, totdep) 

integer channel, startlines(6), endlines(6), skip, 

+ depth(S), curfile, totdep, filenum 

real a, b, mvu, mvl, vl, vu, hvdcout, x, y, sumx, sumy, 

+ t( 1 00), u( 100), v( 100), sumtsq, sumty, sqd, 

+ sqda, sqdb, cov, del, sumxsq 

character* 6 voltrange 
character* 12 filenames(6) 
sumx - 0 
sumy - 0 
sumtsq « 0 
sumty - 0 
sumxsq - 0 
curfile - 1 
totdep - 0 

500 open (1, file-filenames(curfile), form- ’formatted’) 
rewind (1) 

skip - startlines(curfile) - 1 
do 1000 counter - 1, skip 

read (1, *) channel, hvdcout, mvu, mvl, vl, vu 
1000 continue 

depth(curfile) - endlines(curfile) - startlines(curf0e) + 1 
do 1010 counter - (totdep + 1), (totdep + depth(curfile)) 
read (1, *) channel, hvdcout, mvu, mvl, vl, vu 
if (voltrange ,eq. ’lower’) then 
x - mvl 
y » vl 
else 

x - mvu 
y- vu 

endif 

sumx - sumx + x 
sumy - sumy + y 
sumxsq - sumxsq + x*x 
u(int(counter)) - x 
v(int(counter)) - y 
1010 continue 

close (1) 

totdep - totdep + depth(curfUe) 
curfile - curfile + 1 
if (curfile .le. filenum) go to 500 
do 1020 i - 1, totdep 

t(i) - u(i) - sumx / totdep 
sumtsq - sumtsq -t- t(i)*t(i) 
sumty - sumty + t(i)*v(i) 

1020 continue 

a -sumty /sumtsq 
b - (sumy - a*sumx) / totdep 
call sqd2 (a, b, totdep, u, v, sqd) 
del - totdep * sumxsq - (sumx)**2 
sqdb - sqrt (sumxsq / del * sqd * totdep / (totdep - 2)) 
sqda - sqrt (totdep / del * sqd * totdep / (totdep • 2)) 
cov - - sumx / del 
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return 

end 

c sqd2 calculates the square deviation by definition 

c sqd- (1 / n) • sum(Ax + B - y)**2 

subroutine sqd2 (a, b, n, u, v, sqd) 
integer n 

real a, b, u(20), v(20), sqd 
sqd-0 

do 2000 i- 1, n 

sqd - sqd + (a*u(i) + b - v(i))**2 
2000 continue 

sqd - sqd / n 
return 
end 

c The MAIN program body. 

c It combines the files that are corresponding to the same 

c temperature, and produces four lines of output per each 

c combination. The program uses two driver files: 

c mvs-driver.dat and mvs-tc-driver.dat. The lsq2 subroutine 

c has been modified in order to be able to combine several 

c files for analysis. 

integer stll(6), edll(6), stl2(6), edl2(6), pairsl, pairsu, 

+ filnum, i 

real al, bl, sqdl, sqdal, sqdbl, covl, 

+ au, bu, sqdu, sqdau, sqdbu, covu 

character*6 voltrange 
character* 12 filenames(6), mfile, file 
open (2, file- ’mvs-driver.dat’, form- ’formatted’) 
open (3, file»’mvs-tc-out.dat’, form-’formatted’) 
open (4, file- ’mvs-tc-driver.dat’, form- ’formatted’) 

50 i - 0 

read (4, *) mfile 
100 i-i + 1 

read (4, *) filenames(i) 
if (filenames(i) .ne. ’*’) go to 100 
filnum - i - 1 

do 200 i - 1, filnum 
rewind (2) 

300 read (2, *) file, stll(i), edll(i), 

+ stl2(i), edl2(i) 

if (file .ne. filenames(i)) go to 300 
200 continue 

print *, ’Processing master-file: ’, mfile, ’...’ 
voltrange - ’lower’ 

call lsq2 (filenames, filnum, voltrange, stll, edit, 

+ al, bl, sqdl, sqdal, sqdbl, covl, pairsl) 

voltrange - ’upper’ 

call lsq2 (filenames, filnum, voltrange, stll, edll, 

+ au, bu, sqdu, sqdau, sqdbu, covu, pairsu) 

write (3, fmt-2000) mfile, ’chrl’, ’mvl-vl’, al, bl, 

+ sqdl, sqdal, sqdbl, covl, pairsl 

write (3, fmt-2000) mfile, ’chrl’, ’mvu-vu’, au, bu, 

+ sqdu, sqdau, sqdbu, covu, pairsu 

c analyze the second channel range ... 

if (stl2(l) .eq. 0) go to 500 
voltrange - ’lower’ 

call lsq2 (filenames, filnum, voltrange, stl2, edl2, 

+ al, bl, sqdl, sqdal, sqdbl, covl, pairsl) 

voltrange - ’upper’ 
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call lsq2 (filenames, filnum, voltrange, stl2, edl2> 

+ au, bu, sqdu, sqdau, sqdbu, covu, pairsu) 

write (3, fmt-2000) mfile, ’chr2’, ’mvl-vl’, al, bl, 

+ sqdl, sqdal, sqdbl, covl, pairsl 

write (3, fmt-2000) mfile, ’chr2\ ’mvu-vu’, au, bu, 

+ sqdu, sqdau, sqdbu, covu, pairsu 

500 if (mfile .ne. ’f2-e2.dat’) go to 50 
close (2) 
close (3) 

print *, ’File: mvs-tc-out.dat created.’ 

2000 fcamat (al2, 02, a4, tr2, a6, tr2, f9.5, tr2, gl2.6, tr2, 

+ gl 1.6, tr2, gl2.6, tr2, gl2.6, tr2, gl2.6, tr2, i2) 

end 


MEMORANDUM 


Massachusetts Institute of Technology 
Center for Space Research 

September 24, 1987 


VOYAGER MEMORANDUM # 161: 

From: Ralph L. McNutt, 

To: Voyager Internal 

Subject: Modulator Calibrations (Modcal or MVM): Interpretation 


normal Operation 

As part of the routine calibration sequence on the Voyager Plasma Science (PLS) experi- 
ment, there is a digitized readout of the low voltag e signal which is proportional to the poten- 
tial on the modulator grids. This is accompUshetTvia a set of Modulator Voltage Monitors; 
these measurements are referred to by MVM on Summary or EDR tape logs. In many of the 
engineering notebooks they are referred to as "Modcal" measurements or mode (as compared 
to "Curcal" or current calibration measurements. In producing the modulator voltage, the vol- 
tage is stepped through one decade from 60 to 600 volts and then through the same steps with 
a extra gain of x 10 switched into the circuit to sweep from 600 to 6000 volts. To this vol- 
tage is added an offset of 50 volts. Hence, the L and M modes sweep from +10 V to +5950 
V and the E2 mode sweeps from -10 V to -5950 V (the El mode sweeps from -10 V to — 
140 V by using a different step size). Thus, the true voltage ranges across 2.8 decades. i 

The voltage monitors measure a voltage which is a constant fraction of the voltage swept 
which does not include the factor of 10 gain change. The monitors measure a voltage 
between 0.1 and 10 volts which is then input into a fast A/D converter. This voltage is then 
converted to a binary number from which the highest order bit is discarded. Thus, as the low 
voltage sweeps between 0.1 and 1.0 the digital output sweeps from 0 to 255; as the voltage 
sweeps from 1.0 to 10.0 volts, the output again sweeps from 0 to 255, although retention of 
the highest order bit would give a digital output from 256 to 511. Hence, the voltage monitor 
output can correspond to 4 orders of magnitude although, again, the factor of 10 gain change 
is not explicitly incorporated in this measurement To know which decade is actually being 
read out, one must either rely on a "proper sequence" (the case for nominal operation for 
which the MVM readout presumably yields voltages close to the nominal values) or an "edu- 
cated guess" (for interpreting the MVM readouts from the PLS instrument on Voyager 1 after 
the failure). 


nominal voltage Levels 

The nominal conversion between channel number (or "step number") is employed in all 
analysis (both in the VGRANL and MJSANL software packages). The conversion between 
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channel number and voltage is discussed in Bridge et al. [1977]. Let n be the channel 
number, which can take on values of 1 to 16 for the El, E2, and L modes and 1 to 128 for 
the M mode. Let M be an integer with values of 8 for the E2 and L modes, 32 for El and 64 
for M. Then the potential <}>„ at the lower edge ("bottom") of channel n (also referred to~) 
below as V L ^y), which is equal to the voltage at the upper edge ("top") of channel n + 1 j 2 
(also referred to below as V v HV ), (the upper edge potentials follow from considering "chan- ; ■ 
nel" 129 for the M mode and 17 for the other modes), is given by J 

4> n s 60' icy* ^ — 50 (1) 

where <J>„ is in volts. This formula only yields magnitudes, of course, the actual potentials in 
the El and E2 modes being negative. 

It is important to note that the 50 volt offset makes the channel spacing only quasi- 
logarithmic; at lower channel numbers, the channel widths are larger than if they were loga- 
rithmically spaced. 


NOMINAL FAST A/D VOLTAGE CONVERSION 

Information on the analog- to-digital converter is contained in Anton Mavretic’s lab note- 
book labeled book 2, number 254 (starting date of July, 1974). On page 142 (dated 7/25/75), 
we find for analog voltage v running from 0. 1 to 10, the digital count number N running from 
0 to 127. Anton gives the constitutive relation as 

(127 - AO T 

v-V„/ ’ <» 

and values V REF » 10 V, T = » 1/ 57.6 kHz « 17.36 | is, x => 482.55 |xs. Transforming to 

base 10, this yields T~ 

(127 - AQ 

V = 10 x 10 64 00008 

where the denominator in the exponent is usually approximated as the integer 64. 

On page 143, formulas are given for the fast A/D conversion. Anton lists 


( 2 ) 



v = V REF e T , 1.0 < v < 10. 

(3a) 

and 

(255 - M T 

v = V REF /\0e T , 0.1 < v < 1.0. 

(3b) 

with T=fi£, ad = 1/ 230 - 4 ^ = 4 34 x = 482 55 V*- Transforming to base 

10, this 

yields the combination 

(255 - AO 

v = 10 x 10 256 , 1.0 < v < 10. 

(4a) 

and 

(255 -AO 

v = 10 256 , 0.1 < v < 1.0. 

(4b) 


Although not explicitly stated, these are apparently the nominal, rather than measured, 

- 2 - 



conversion formulas. 


NOMINAL MONITOR OPERATION 

In the usual MVM mode, the data numbers (DN) produced by running the low voltage 
monitor outputs through the fast A/H converter are returned as a function of channel number. 
In the L and M modes, the DNs corresponding t o V^ a re returned in the space occupied by A 
cup plasma measurements in the usual plasma measurement mode; similarly, the DNs 
corresponding to V H are returned in the space occupied by the B cup plasma measurements. 
The C cup and D cup contain currents corresponding to normal plasma currents but with DN 
values obtained from th e, fast A /D converter. Hence, for the L and M modes, there is a meas- 
ure of the contiguity of neighboring channels (gaps and/or overlaps in coverage). 

For the El and E2 modes, there are only 16 slots available in the data sequence for each 
mode, so compromises were made in deciding which DN values would be included in the 
telemetry stream. Fo r odd c hannel numbers (counting from 1 to 16), the DN corresponding to 
Vr. is returned. For even channel numbers, the DN corresponding to V v is telemetered. 
Hence, both absolute value and contiguity are checked, but only for every - jther channel. For 
example, there are MVM measurements for V Ll , V w , V u , V UA , V^, etc., but there is no cali- 
bration information for the boundaries between channels 1 and 2, 3 and 4, etc. 

In nominal operation, note that we should have Vu. , = v L . Indication of where the 
V(js and V^s come in the data is given in the VOYPRT appended to Voyager Memorandum 
158 as Figure 6. 


RELATIONSHIP BETWEEN MODULATOR POTENTIAL 
AND MONITOR VOLTAGES 

Figure 1 (from J. Binsack) shows a simplified block diagram of the high voltage modulator 
circuit for the PLS instrument; indicated on the figure is the location of the suspected failed 
circuit in the Voyager 1 PLS instrument The low voltage monitor voltages are measured at 
the point indicated with the arrow and labeled = CAL. This low voltage is (nominally) 
directly proportional to the high voltage applied to the modulator grids (at least in the linear 
regime of the amplifiers). The DC gain of the Cockcroft-Walton amplifiers ("C. W.") and the 
AC gain of the transformer in the feedback loop are unknown as of this writing, so a theoreti- 
cal value of the proportionality factor is not known (but could be derived from the appropriate 
circuit diagrams, J. Binsack, private communication, 1987). Data shown in Appendix A (see 
below) suggests that this factor is - 60. 

The subcontractor responsible for the high voltage power supply, Matrix Research and 
Development Corporation, measured both the voltage to ground from the modulator grid leads 
and the low voltage monitor outputs for all three modulators (on SN001, SN002, and SN003) 
for all channels in all modes. * Work has begun to combine this information with other infor- 
mation to develop an algorithm to link the measured DNs in the MVM mode with the actual 
modulator potential. It is important to note that all modulator potential measurements were 
actually made with high impedance probes to the modulator grid leads, terminated with capa- 
citances to approximate the electrical characteristics of the grids themselves (R. Butler, private 
communication, 1987). In no case were potential measurements made on the grids themselves 
with the grids in place in the cups; this is not possible as a probe cannot physically be 
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inserted into the assembly. 

Before the PLS instruments were to delivered to JPL for integration with the Voyager 
spacecraft, extensive calibration tests were run at the MIT Laboratory for Space Experiments 
using the Bench Checkout Equipment (BCE). The modulator was stepped through all chan- 
nels in all modes and the DNs (or "decimal counts") were recorded as a function of step 
number, equal to the channel number minus one (channel 1 corresponds to step 0). 

In Appendix A, the reference counts are listed as a function of step number for each mode 
(labeled for the prototype rather than for one of the flight units). Corresponding voltages 
from the voltage monitors are also given. Note that these are the counts measured with the 
BCE^ For example, the nominal value of V L yy for step 0 is +10 V in L and M modes and 
-10 V in El and E2 modes. The reference counts for the L, El, E2, and M modes are, y 
respectively, 54, 48, 50, 48. The fraction values for V L are, presumably, the measured low 
voltage monitor outputs (they are not nominal values as the nominal value should be the same 
in for all of these). The applicable equation relating the monitor voltages and the reference 
counts is (4b). However, applying (4b) to values of 0.172, 0.163, 0.166, 0.163, yields DNs of 
59.3, 53.3, 55.3, and 53.3, respectively. This suggests that the theoretical A/D algorithms are 
off (if the reference counts and monitor voltages were really measured simultaneously). Note 
that (4b) is off by about a factor of lO 5 ' 256 = 1.046, but not quite. 

The agreement is better at the other end of the scale for Vy of step 15 of E2 and L and step 
127 of M, the counts are 248 and the monitor voltage is 9.314 V. Using 9.314 volts and 
equation (4a) yields a DN of 247.1. Again^if th e counts a nd monitor voltages were measured 
at the same time (which is not known), there is some nonlineari ty in the A/D converter, not 
"currently accounted for. A spot check suggests that a simple DC offset voltage also does not 
account for the discrepancy. 

~ I have searched through all of Anton’s notebooks stored in Herb’s office and through all of 
the files in N-52 and have not been able to locate any documentation (other than that men- 
tioned above) on the performance of the A/D converter. To put all of this in perspective, at 
the higher voltages in the M mode, the reference counts for Vy and V L are separated by 4 
DN. This implies that a determination of the true potentials on the modulator grids to an 
accuracy of better than 1% should, in principal, be possible. 

Attached as Appendix B are similar modulator calibration data from the Flight 1 unit (on 
Voyager 1, i.e., SN002). Most of the measurements are dated 4/22/77, less than 5 months 
before the launch. Unfortunately, simultaneous high voltage, monitor voltage, and DN were 
not recorded. It is worth noting that for step 127 of the M mode, the measured value of Vu 
was 6010 volts, as compared to a nominal value of 5950 volts, a deviation from nominal of 
1%. This suggests that such deviations are not uncommon. 

MODULATOR voltages for the VOYAGER 1 PLS 
instrument after the Failure 

The assumed circuit failure on the Voyager 1 PLS instrument pegged the value of V L at 
some value greater than 6000 V, and caused the values of Vy to increase with step number, 
but in a spurious way (J. Binsack, private communication, 1987), i.e., the phase of the modu- 
lated current is shifted 180°. Hence, the solar wind can no longer be detected in the normal 
plasma mode because the current is detected synchronously, and the failure introduced the 
incorrect phasing. 
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If the value of V L is pegged to a large value, then the solar wind is always excluded from 
the detector for half of a measurement. If the value of Vy is low enough, then the solar wind 
should be admitted for the other half, unless Vy steps to a large enough value to also exclude 
the solar wind. 

The DC return measurements from the Voyager 1 instrument after the failure show changes 
in all 16 channels which vary with time; this presumably gives some measure of the solar 
wind flux (see Voyager Memorandum 157). There is a "break" in the DN numbers in the DC 
return mode above channel 8, suggesting that Vy for the first 8 channels (in the L mode) is 
always low enough to allow the solar wind access to the collector plate. The DNs decrease 
for the top 8 channels suggesting that more of the solar wind is gradually being excluded, but 
the lack of variation in this pattern could mean that this is another instrumental effect instead. 
The interpretation is also complicated by the fact that there is a factor of 10 gain change inter- 
nally in going from the voltages for channel 8 to those of channel 9, so the break at this point 
may actually be related to this change. 

From the L mode spectrum in Figure 6 of Voyager Memorandum 158 (referred to previ- 
ously), the MVM data numbers change from 190 to 93 for the first 8 channels of the L mode 
(at 1980-346/0103:23.735) with a "fold" occurring between channels 2 and 3. If the DNs can 
be interpreted as before the failure, then the fold could occur at -6, 60, 600, or 6000 V. The 
value of 6000 V must be excluded on the basis of the detection of the solar wind. At channel 
8, the DN of 93 could be -14 or 140 V; higher values are again excluded because the solar 
wind is detected. 

In the top 8 channels, the DNs run from 90 to 136. If no "folds" are present (and none are 
apparent), this could correspond to -140 V to 210 V or some multiple of 10 times these lim- 
its. These potentials would allow the solar wind to be admitted in all channels, but for speeds 
of 400 km/s, similar to those detected by Voyager 2 at these distances, no decrease in DC 
current with channel number should be present. Multiplying the limits by a factor of 10 
would exclude the solar wind from the top channel, unless its speed exceeded -600 km/s, and 
speeds of this magnitude were not detected just prior to the failure, so something is not con- 
sistent. 

To understand the post-failure measurements by the instrument on Voyager 1, a more 
detailed investigation is clearly required. Such an investigation is now underway. 
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